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ABSTRACT 

Silicate dnst grains in the interstellar medinm are known to be mostly amor- 
phons, yet crystalline silicate grains have been observed in many long-period 
comets and in protoplanetary disks. Annealing of amorphons silicate grains into 
crystalline grains reqnires temperatnres > 1000 K, bnt exposnre of dnst grains in 
comets to snch high temperatnres is incompatible with the generally low temper¬ 
atnres experienced by comets. This has led to the proposal of models in which 
dnst grains were thermally processed near the protoSnn, then nnderwent consid¬ 
erable radial transport nntil they reached the gas giant planet region where the 
long-period comets originated. We hypothesize instead that silicate dnst grains 
were annealed in situ, by shock waves triggered by gravitational instabilities. We 
assnme a shock speed of 5kms“^, a plansible valne for shocks driven by gravita¬ 
tional instabilities. We calcnlate the peak temperatnres of micron and snbmicron 
amorphons pyroxene grains of chondritic composition nnder conditions typical in 
protoplanetary disks at 5 - 10 AU. Onr results also apply to chondritic amor¬ 
phous olivine grains. We show that in situ thermal annealing of submicron and 
micron-sized silicate dust grains can occur, obviating the need for large-scale 
radial transport. 

Subject headings: comets: general — dust — solar system: formation — shock 
waves 
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1. Introduction 


Comets are collections of the dust grains and ices that resided in the solar nebula. 
The dust grains and ices in comets do not appear to have been subject to any processing 
mechanisms since the comets formed at the birth of our solar system (Hanner 1999; Brucato 
et ah 1999). Cometary material therefore is thought to contain pristine solar nebula material 
and to represent the physical conditions of the solar nebula at the birth of the comet, 
including temperatures <50K (Meier & Owen 1999). Complicating this interpretation is 
the spectroscopic evidence of silicate grains processed at temperatures > 1000 K existing in 
comets. Crystalline silicates with distinct mid- and far-IR spectral resonance features can 
form a signihcant component of dust (30 - 50%) in comet comae (Wooden et ah 1999; Barker 
et ah 2002). The majority of the crystalline silicate grains seen in many comets must have 
been condensed or annealed (presumably by thermal processing) in the solar nebula. As 
pointed out by Irvine et ah (2000) and Wooden et ah (1999), despite the observations of 
crystalline olivines [(Mg 3 ,Fea;_i) 2 Si 04 ] and pyroxenes [(Mg^Fea,_i)Si 03 )] {x ranges between 1 
[Mg-pure] and 0 [Fe-Pure]) in circumstellar environments (Waters & Waelkens 1998; Malfait 
et ah 1998), crystalline silicates are essentially unseen in the interstellar medium (< 5%; Li 
& Draine 2001). However, it is not clear how silicate grains could be heated to ~ 1000 K in 
the cold, comet-forming environment (~ 5 — 100 AU). 

In the solar nebula, the high temperatures > 1000 K needed to anneal amorphous in¬ 
terstellar silicate grains, or to condense crystalline silicate grains directly from the gas phase 
(> 1300 K) were achieved in steady-state only at distances < 1 AU from the protoSun at 
early times (e.g.. Bell et ah 1997). Proposed mechanisms for transporting crystalline sili¬ 
cates to the comet forming zones include radial diffusion from vigorous thermal convection 
(Nuth 2001; Hill et al. 2001), or outflows driven by reconnecting magnetic held lines (“X- 
winds” [Shu et al. 1996]). Although X-winds have been invoked as the cause of chondrule 
formation (mm-sized spheres of silicate rock found in meteorites [see Jones et al. 2000]), 
there is compelling evidence (chondrule rim thickness, chondrule-matrix complementarity) 
that chondrules formed in situ, in the asteroid belt region (see summary in Desch & Connolly 
2002; hereafter DC02), obviating the need for radial transport. The similarities between the 
annealing of silicate grains and the melting of chondrules have been pointed out by Grady et 
al. (2000) and Nuth (2001). These similarities between chondrule formation and silicate an¬ 
nealing lead us to hypothesize that interstellar, amorphous silicate grains were also annealed 
in situ, at 5 - 10 AU, by shock waves in the solar nebula. 

In this Letter, we calculate the thermal histories of submicron silicate grains passing 
through nebular shocks, using the code of DC02 that was used to calculate the thermal 
histories of chondrules. We conclude that silicate grains were likely annealed in situ in the 
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gas giant planet region by nebular shocks, and we attribute these shocks to disk gravitational 
instabilities. 


2. Annealing in Shocks 
2.1. Silicate Grains and their Emissivities 

Dust grains of radius Op are heated in shocks by thermal exchange with the gas and by 
frictional gas-drag heating, and cool at a rate iTTOp Qp cT^, where a is the Stefan-Boltzmann 
constant, and Tp and Qp are the particle’s temperature and emissivity. We calculate the 
emissivities of amorphous pyroxene grains of a chondritic composition [i.e., (Mgp 5 Feo. 5 )Si 03 ] 
ranging from 0.1 /rm to 1.0 /rm in radius, using optical constants from Dorschner et al. (1995). 
Our chosen grain size and composition correspond to grain properties determined from the 
study of IDPs (Bradley et al. 1999; Bradley et al. 1996), solar system comets (Barker et 
al. 2002; Mason et al. 2001; Schulze, Kissel, & Jessberger 1997), and circumstellar disks 
(Bouwman et al. 2000; Malfait et al. 1998). The opacities of amorphous pyroxene and 
amorphous olivine are similar enough (< 5% discrepancies) that our results can be extended 
to olivine as well. The wavelength-dependent emissivity is found as a function of particle 
size using Mie theory. 

The Silicate Evolution Index, or SEI (§2.3), which describes the degree of annealing 
of silicate material, strictly refers only to Mg-rich silicates, not the chondritic composition 
chosen here. ISO SWS found Mg-rich {x > 0.9) crystalline silicates in comet Hale-Bopp 
and in the circumstellar environments around Herbig Ae/Be stars. In the experiments of 
Hallenbeck et al. (2000), only Mg-rich crystals (forsterite, enstatite) were annealed. We 
chose an Fe-rich, chondritic compostion silicate to generate conservative estimates of the 
peak temperature of grains. Grains heated in shocks by frictional drag heating and thermal 
exchange with the gas cool by emitting thermal radiation at ~ 3 /rm. The opacity at 3 /rm 
is much higher for grains of chondritic composition than for pure Mg-rich compositions. 
Accordingly, Mg-rich grains cool less efficiently and reach even higher temperatures than the 
grains considered here. Mg-rich grains would be more easily annealed than the grains we 
consider. 


2.2. Calculation of Temperatures 

To calculate the temperatures of particles overtaken by nebular shocks, we use the code 
of DC02, to which the reader is referred for all computational details. As described in DC02, 
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shocks waves heat, compress, and accelerate gas almost instantaneously. Particles in the gas 
come into dynamical equilibrium with the gas through frictional forces, which can heat the 
particles. Particles are also kept hot by thermal exchange with the hot gas, and by absorbing 
the radiation emitted by other particles. Absorption of thermal radiation heats the particles 
even before they reach the shock front. For the purpose of calculating the radiation held, 
we assume a dust opacity of 4cm^ per gram of gas, due to small (< 0.1/xm) particles in 
thermal equilibrium with the gas. Particles cool by emitting radiation according to their 
emissivities (§2.1). In practice, the thermal histories of particles and gas are calculated 
assuming a best guess for the radiation held. An updated radiation held is then calculated, 
and the thermal histories recalculated. Iterations are continued until the temperatures of 
particles at each location are converged to within < 10 K, which takes on the order of 10^ 
iterations. Since particle temperatures only increase with each iteration, the convergence 
of the solution is well-dehned, and the calculated SEI are lower limits. All details of the 
calculation are identical to DC02, with the exception that gas properties appropriate to 
the 5-10 AU region are used (§2.4), and the particle radius is ~ 1/im. The thermal 
and dynamical timescales scale with particle radius, so this small particle radius makes the 
calculation a factor of 10^ more computationally intensive than the chondrule program. Each 
calculation presented here took a week on a PC with dual 833 MHz processors, so the SEI 
is calculated for a small number of parameter combinations. 


2.3. Annealing 

The annealing of amorphous silicate grains into crystalline grains has been studied ex¬ 
perimentally by Hallenbeck et ah (2000), who have developed a predictive quantity called the 
Silicate Evolution Index (SEI). This quantity is essentially an integration of the temperature 
T over the thermal history of the silicate particle, and is meant to predict the change in the 
spectra of silicate particles as they anneal. At temperatures < 1067 K, the spectra evolve in 
a two-step process: as grains spend a time t at high temperatures, their spectra evolve, then 
stall, then continue evolving. Eor temperatures that exceed 1067 K, there is no stall in the 
evolution of the spectra; in this case, SEI si 1 -|- 5.88 x 10®® t exp (—177,430/T), with t in 
seconds. Grains with SEI > 2 are evolving toward crystalline, and SEI > 100 are considered 
completely crystalline. Using this formula, we can estimate the peak temperatures grains 
must reach to be annealed in the shock. In our parameter studies, particles are found to be 
within 20 K of their peak temperatures for < 0.1s (see § 2.4). Values of the SEI exceeding 
10 (moderately annealed) therefore require peak temperatures ^ 1150 K. Empirically (Table 
1), annealing was found to require peak temperatures ~ 1200 K. 
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2.4. Results 

The temperature of a 1.0 pm grain during the ten minutes before and ten minutes after 
being overrun by a 5kms“^ shock is shown in Fig. la. As discussed in §3.1, we consider 
this a likely speed of shocks at 5 - 10 AU due to gravitational instabilities. The ambient gas 
density in this example is Pg = 5 x 10“^*^ g cm“^. Most of the annealing of this pyroxene grain 
will take place while the particle is at its peak temperature, 1200 K, at which it remains for 
~ 0.1s, as shown in Fig. 16. Integrating over the thermal history of this particle yields a 
value of the SEI ^ 8, meaning that this particle is marginally annealed by passage of a shock 
through gas of this density. 

The peak temperatures of initially amorphous pyroxene grains overrun by 5kms“^ 
shocks, and the corresponding values of the SEI, are listed in Table 1 for other combinations 
of gas density and particle size. The particle radii, 0.1,0.25,0.5 and 1pm, were chosen to 
match the inferred sizes of crystalline silicate grains in comets. The ambient gas densities, 
0.3,1,2, 3, 5 and 10 x lO^^'^gcm”^, were chosen to correspond to the densities computed in 
the a disk models of Bell et ah (1997), parameterized by the mass accretion rate M through 
the disk and by the standard Shakura-Sunyaev turbulence parameter a. For example, if 
M = 10“® Mq yr“^ (the median mass accretion rate in T Tauri systems; Gullbring et al. 
1998), and a = 10“^, then pg = 8.4 x 10“^°gcm“^ at 5 AU, and 3.0 x 10“^°gcm“^ at 10 
AU. For a = 10“^, the values are 2.8 x 10“^^gcm“^ at 5 AU and 1.2 x 10“^^gcm“^ at 10 
AU. The post-shock gas densities are greater than the ambient density by a constant factor, 
~ 6, and the post-shock temperature are set largely by the shock speed only. For many of 
the cases seen in Table 1, we hnd SEI >1. Therefore, shock-annealing of silicate grains in 
the 5-10 AU region is evidently possible. 


3. Discussion 

Annealing of silicate grains can be accomplished by shocks in the outer solar nebula (at 
10 AU). Smaller grains are more easily annealed due to their poorer radiative efficiencies. 
These results are not affected by the presence of ice mantles, because such mantles would 
sublimate before reaching the shock front. Water ice is at least as good an absorber as 
silicates of the thermal radiation emanating from the shock front, and dust grains can easily 
absorb enough energy from the radiation held to heat above 150 K. Figure la shows that 
a dust grain attains temperatures of > 1000 K for > 10 minutes. Since the grain is in 
radiative equilibrium, the grain is absorbing and emitting thermal radiation at a rate of 
hTTUp Qp uTp. If the grain is assumed to be made entirely of water ice with latent heat of 
sublimation 2.83 x 10^'^ erg g“^ and radiative absorption efficiency Q = 0.1, it can easily 
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be shown that the ice will completely evaporate in < 0.2 s. This evaporation time is much 
less than the 10 minutes the particle stays at temperatures > 1000 K. Shocks raise grains to 
high temperatures whether or not they have ice mantles. Therefore, provided shock speeds 
~ 5kms“^ can be attained, annealing of silicates in situ is very likely. 

Shocks are very likely to occur early in the evolution of the solar nebula. The accre¬ 
tion disk models of Bell et al. (1997) predict that T Tauri disks with mass accretion rates 
lO“®M 0 yr“^ and a = 10“^ are locally gravitationally unstable beyond 11 AU, and disks 
with a = 10“"^ are locally unstable beyond 6 AU. Since mass accretion rates ~ 10“® Mq yr“^ 
are likely to persist for many Myr in T Tauri systems (Gullbring et al. 1998), gravitational 
instabilities are likely to persist in the outer regions of protoplanetary disks. Since the gravi¬ 
tational instabilities manifest themselves as global modes (Laughlin & Rozyczka 1996), they 
are likely to affect the bulk of the disk, even where it is locally stable. 

How the gravitational instabilities in a marginally unstable disk manifest themselves 
has been explored numerically by Boss (2000, 2001a,b), who discovered the development 
of gravitationally bound clumps and highly non-axisymmetric density structures associated 
with the clumps in such disks. These structures would be overdense by orders of magnitude, 
oriented more or less radially, extending inward from the clumps and co-rotating with them 
(see Figure 1 of Boss [2000]). Gas orbiting closer to the Sun would orbit faster than these 
density structures and collide more or less normally with them, at relative velocities com¬ 
parable to the difference in Keplerian angular velocities between the orbiting gas and the 
clump driving the structure. For example, a clump at 10 AU would drive a density structure 
with an orbital period ~ 30 yr; material orbiting every 15 yr at 6 AU would slam into this 
structure at a relative velocity ~ 6 kms“^. 

If nebular shocks existed in the outer solar nebula, whatever their source, we hud that 
submicron silicate grains could readily be annealed by 5kms“^ shocks, if the gas densities 
were ^ 10 “^°gcm“^. These densities are plausible for the 5-10 AU region of the solar 
nebula, according to the a disk models of Bell et al. (1997), but not beyond 10 AU. Therefore, 
comets forming in the Kuiper Belt Region (greater than 35 AU) should not contain crystalline 
silicate dust annealed by shocks. Radial diffusion by turbulence is unlikely to change this 
result. 

Within the time a comet takes to form, t ~ 10^ yr (Weidenschilling 1997), turbulent 
diffusion mixes gas across distances ~ where C is the sound speed and H the 

scale height of the nebula. For a ~ 10“^, the mixing distance over 10® years at 10 AU is 
<5AU; over 10® years it is <15AU. 

The vast majority of short-period comets, which come from the Edgeworth-Kuiper belt 
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and which mostly formed beyond 35 AU, should contain amorphous grains only. The few 
short-period comets with crystalline silicates must have been scattered to the Edgeworth- 
Kuiper belt after forming at 5 - 10 AU, which is possible dynamically (Duncan & Levison 
1997). However, most comets formed at 5 - 10 AU will instead be scattered into the Oort 
cloud and be observed entering the solar system today as long-period comets. It is in long- 
period comets that crystallinity should be expected. 

The observational data conhrm the prediction that crystalline silicates should be found 
almost exclusively in long-period comets. Fits to Hale-Bopp’s (C/1995 01) 11.2 /rm features 
have led to the determinations that, even pre-perihelion, si 20% of all silicates (Hanner et al. 
1997, Wooden et al. 1999), and ~ 30% of the olivine silicates (Harker et al. 2002) must be in 
crystalline form. According to Yanamadra-Fisher & Hanner (1999), comets with comparable 
peaks at 10 and 11.2 /rm, such as P/Halley, Bradheld 1987 XXIX, Mueller, Levy 1990 
XX, and Hale-Bopp, are well fit by a combination of 20% crystalline and 80% amorphous 
silicates. The crystalline silicate feature at 11.2/rm has also been observed in Hyakutake 
(C/1996 B2; Mason et al. 1998), and McNaught-Hartley (C/1999 Tl; Lynch, Russell & Sitko 
2001). Not all long-period comets have crystalline silicates —Hanner et al. (1994) found no 
evidence for annealing in Kohoutek 1973 XH, Austin 1990 V, Okazaki-Levy-Rudenko 1989 
XIX, and Wilson 1987 VH —but most do. On the other hand, only in one short-period comet, 
103/P Hartley 2, have crystalline silicates been observed (Crovisier et al. 1999). Data are 
inconclusive or negative in others (Harker et al. 1999; Hanner et al. 1996), including comets 
P/Brorsen-Metcalf 1989 X (Hanner et al. 1994), P/Churyumov-Gerasimenko (Hanner et 
al. 1985), P/Grigg-Skjellerup (Hanner et al. 1984), and P/Schaumasse (Hanner et al. 1996), 
each indistinguishable from blackbodies (Hanner et al. 1994). Winds hypothetically carrying 
annealed grains from the inner solar system would preferentially seed Oort-cloud comets with 
crystalline silicates, but these winds could possibly also deposit crystalline silicates in Kuiper- 
belt comets as well. The low gas densities and orbital velocities of the shock model rule out 
this possibility for comets that form beyond ~ 10 AU. 

If shocks annealed silicate grains in situ, the possibility exists that the nebula gas could 
be chemically processed by the shock as well. For example, Kress et al. (2001) suggest that 
nebula shocks could produce nitriles that are incorporated into nitrogen-bearing organics in 
meteorites. Observed correlations of shock-produced, disequilibrium chemicals like nitriles, 
with silicate crystallinity, would strengthen the case for shock annealing. 

Comets, nearly all of them long-period comets, despite forming at very cold temper¬ 
atures, contain crystalline silicate grains that imply temperatures ^llOOK, if only for a 
fraction of a second. These conditions are met in nebular shocks. Gravitational instabilities 
make shocks a likely event in the hrst few Myr of the solar nebula, but whatever the source. 



silicate grains can be annealed by shocks in the 5-10 AU region of the solar nebula. Shock¬ 
annealing obviates the need for large-scale radial transport of annealed grains from the inner 
solar nebula to the 5-10 AU region where the Oort-cloud comets formed. 
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Table 1. Peak Temperature and SEI of Annealed Grains 


p (g cm 

0.1 /im 

grain 
0.25 pm. 

radius 

0.5 pm 

1.0 pm 

0.3 X 10-^° 

1100 K 
4.0 X 10-° 



670 K 

4.9 X 10-° 

1.0 X 10-^° 

1305 K 

1095 K 

990 K 

920 K 


6.6 X 10^ 

2.2 X 10-° 

3.6 X 10-° 

1.1 X 10-® 

2.0 X 10-1° 

1405 K 

1210 K 

1110 K 

1035 K 


3.3 X 10° 

4.1 X 10° 

7.3 X 10-° 

1.1 X 10-® 

3.0 X 10-1° 



1175 K 
4.3 X 10-2 

1110 K 

2.7 X 10-° 

5.0 X 10-1° 



1255 K 

3.8 X 102 

1200 K 

8.3 X 10° 

10.0 X 10-1° 

1650 K 
2.4 X 10i° 



1300 K 

6.0 X 10° 
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Fig. 1.— Temperatures of a dust particle (solid line) and gas (dashed line) in a nebula shock. 
The ambient gas density, 5 x 10“^'^gcm“^, and shock speed, 5kms“^, are just sufficient to 
anneal this 1.0 /rm-radius grain. The temperatures seen by the particle ten minutes before to 
ten minutes after passing through the shock front are shown in (a). Absorption of thermal 
radiation emitted by warm dust heats the particle in both the pre-shock and post-shock 
regions, and thermal exchange with hot gas heats the particle in the post-shock region. 
A short-lived spike is evident in (b), which displays the temperatures in the two seconds 
surrounding passage through the shock. Friction with the gas adds to the heating during the 
time it takes the particle to come into dynamical equilibrium with the gas, and this extra 
heating to 1200 K leads to annealing of the particle, with SET si 8. 









